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I.  INTRODUCTION 


Two  aspects  of  digital  communication  were  investigated.  Chapters  II,  III  and  IV 
present  an  analysis  determining  the  advantages  of  non-uniform  windowing  in  reducing 
the  effects  of  frequency  offsets  caused  by  doppler  shifts  and  other  sources  in  a  fading 
channel.  Chapters  V,  VI  and  VII  develop  a  simulation  of  adaptive  automatic  repeat 
request  (ARQ)  protocols  to  improve  throughput  performance  through  various  bit 
error  rates. 

A.  EFFECTS  OF  NON-UNIFORM  WINDOWING 

Digital  communication  involves  the  transmission  of  information  in  bit  format. 
M-ary  frequency-shift-keying  (MFSK)  uses  M  =  2^  different  frequencies  to  represent 
M  different  symbols  (each  containing  k  bits  of  information).  A  proper  determination 
of  the  frequency  sent  enables  the  receiver  to  determine  which  symbol  was  transmitted. 

With  the  advent  of  real-time  Fast  Fourier  Transform  (FFT)  processors,  an 
MFSK  receiver  can  be  easily  implemented.  Selected  output  bins  of  the  FFT  will 
correspond  to  the  M  possible  frequencies  provided  the  sampling  rate  and  the  signal 
frequencies  are  closely  related.  In  the  case  developed  here,  a  bin  separation  of  two 
will  be  used  (proper  determination  of  sampling  rate  and  corresponding  bin  locations 
will  be  discussed  in  Chapter  II).  By  comparing  the  magnitudes  of  the  FFT  output 
bins  and  choosing  the  bin  with  the  largest  magnitude,  the  transmitted  symbol  ran  be 
identified. 

Noise  in  the  channel  causes  a  variance  in  the  magnitudes  of  the  FFT  output 
and  an  error  can  occur  when  the  magnitude  of  one  output  bin  is  greater  than  the  bin 
corresponding  to  the  frequency  of  the  signal  sent.  In  the  presence  of  frequency  offsets 


(caused  by  doppler  shifts  or  receiver  oscillator  drift,  for  example),  the  probability 
of  error  is  further  increased  by  a  shifting  of  signal  power  into  other  frequency  bins. 
Frequency  offsets  also  contribute  indirectly  to  erroneous  frequency  components  by 
causing  leakage.  Leakage  is  encountered  in  FFT  analysis  and  is  a  result  of  a  periodic 
signal  not  being  sampled  over  an  integer  number  of  cycles.  A  discontinuity  will  occur 
at  the  endpoints  of  the  signal  and  will  contribute  erroneous  frequency  components  to 
the  FFT.  These  contributions  can  be  seen  in  Figure  1.1(a)  where  the  magnitudes  of 
a  FFT  are  given  for  a  noiseless  signal  experiencing  a  frequency  offset. 


Figure  1.1:  Effects  of  frequency  offset  with  (a)  rectangular,  and  (b)  Hamming 
windowing.  Dynamic  range  is  improved  using  non-uniform  windowing  with  a  loss  of 
signal  detectability. 

Non-uniform  windowing  can  be  used  on  the  sampled  signal  in  order  to  smooth 
out  discontinuities  at  the  endpoints  of  the  sampling  interval  and  help  reduce  the  effects 
of  leakage.  The  reduced  effects  of  leakage  are  balanced  against  signal  attenuation 
caused  by  non-uniform  windowing.  With  less  energy,  signal  detectability  deteriorates 
and  the  noise  variance  has  a  greater  effect  on  proper  signal  decoding.  This  is  seen  in 
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Figure  1.1(b)  where  a  Hamming  window  is  applied  to  the  sampled  signal.  Although 
the  dynamic  range  between  the  two  components  of  interest  has  increased,  signal 
strength  has  decreased.  As  the  frequency  offset  becomes  more  dramatic,  the  benefits 
of  reducing  the  amount  of  leakage  by  non-uniform  windowing  outweighs  the  loss  in 
signal  detectability. 

This  thesis  investigates  the  amount  of  frequency  offset  necessary  to  overcome 
the  performance  loss  due  to  signal  attentuation  in  a  Rician-fading  channel.  The 
results  of  [1]  are  applied  to  a  channel  experiencing  slow  Rician  fading  to  determine 
the  effects  of  fading  on  the  performance  trade-off  between  non-uniform  and  uniform 
windowing.  A  statistical  analysis  is  carried  out  to  determine  the  probability  of  bit 
error  as  a  function  of  sigiial-to-noise  ratios  as  well  as  direct-to-fading  ratios  associated 
with  Rician  fading.  From  these  results,  the  frequency  offset  at  which  non-uniform 
windowing  out-performs  uniform  windowing  is  determined.  It  is  shown  that  as  the 
channel  approaches  Rayleigh  fading  (the  direct-to-fading  ratio  decreases),  the  amount 
of  frequency  offset  necessary  to  justify  the  use  of  non-uniform  windowing  decreases 
despite  its  signal  attenuation  characteristic. 

B.  ADAPTIVE  ARQ 

Having  an  understanding  of  the  expected  bit  error  rate  experienced  by  the 
channel.  Chapters  V  through  VII  investigate  the  possibility  of  improving  throughput 
performance  by  implementing  an  adaptive  automatic  repeat  request  (.ARQ)  scheme 
in  computer  communication. 

Communication  betw'een  computers  can  be  described  using  the  seven  layer  open 
systems  interconnection  (OSI)  model  given  in  Figure  1.2.  Information  from  each  node 
is  first  processed  down  to  the  network  layer  and  into  the  data  link  control  (DLC). 
The  DLC  then  prepares  the  packet  of  information  for  transmission.  Overhead  bits 
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are  appended  to  the  packet  which  include  coding  information  for  error  detection  and 
synchronization.  In  addition,  a  header  and  a  tail  is  attached  to  the  packet  so  that  it 
may  be  identified  by  the  proper  receiving  node.  The  processed  packet,  now  a  frame, 
is  released  from  the  DLC  to  the  physical  interface  for  transmission  over  the  physical 
link.  The  transmitted  frame  is  then  received  at  the  receiving  node’s  physical  interface 
and  a  reverse  process  occurs. 


Figure  1.2:  In  an  OSI  network  communication  representation,  data  packets  are 
passed  from  the  network  layer  to  the  data  link  control  (DLC)  layer  where  they  are 
formatted  and  processed  for  transmission  over  the  physical  link.  ARQ  protocols  are 
carried  out  in  the  DLC  layer. 

As  noted  in  the  previous  section,  noise  on  the  physical  link  sometimes  corrupts 
the  transmitted  frame  and  an  error-detection  scheme  is  implemented  by  the  DLC  to 
recognize  transmission  errors.  When  an  error  occurs,  the  receiving  DLC  must  issue  a 
command  to  the  sender  of  the  frame  to  request  retransmission  of  the  frame. 
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Several  different  ARQ  protocols  currently  exist  and  three  will  be  analyzed  here — 
stop- and- wait,  go-back-N  and  selective- repeat.  An  analysis  of  the  throughput  of  each 
will  be  developed  to  compare  the  performance  of  each  protocol.  From  this  analysis, 
it  is  found  that  different  fixed  frame  lengths  provide  better  performance  for  various 
bit  error  rates.  Although  one  frame  length  may  be  efficient  at  a  given  bit  error  rate, 
that  same  length  may  not  be  as  efficient  if  the  bit  error  rate  changes.  This  is  reflected 
in  Figure  1.3  where,  although  a  frame  length  of  1024  bits  is  most  efficient  at  low 
bit  error  rates,  at  high  bit  error  rates  a  substantial  improvement  in  throughput  is 
achieved  using  a  frame  length  of  128  bits. 


10“  '  ^  "lb'"  "  lb- 

Probability  of  Bit  Error 


Figure  1.3:  Dotted  curves  represent  throughput  of  a  fixed  frame  length  .^RQ  proto¬ 
col  for  various  frame  lengths  L.  At  low  bit  error  rates,  increased  percentage  overhead 
in  small  frames  decreases  throughput.  .At  higher  bit  error  rates,  longer  frames  have  a 
higher  probability  of  error.  The  solid  curve  represents  the  desired  performance  from 
an  adaptive  .ARQ  scheme. 

To  help  improve  performance  through  the  entire  range  of  bit  error  rates,  an 
adaptive  .ARQ  protocol  can  be  implemented.  With  adaptive  frame  lengths,  long 


frames  are  used  at  lower  bit  error  rates  and  shorter  frames  at  higher  bit  error  rates. 
By  adapting  to  the  appropriate  frame  length,  the  throughput  achieved  can  approach 
the  maximum  at  any  BER  of  the  channel  (the  solid  curve  in  Figure  1.3).  The  concept 
of  an  adaptive  ARQ  strategy  is  proposed  in  [2].  This  thesis  presents  the  results  from 
simulations  where  the  length  of  a  frame  can  adapt  to  the  experienced  bit  error  rate 
of  the  channel. 
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II.  RECEIVER  ANALYSIS 


An  analysis  similar  to  that  presented  in  [1]  is  carried  out  where,  in  an  M-ary  FSK 
communications  scheme,  one  of  M  signals  is  sent  during  the  time  interval  0  <  t  <  T. 
Customarily,  M  =  2*^  so  that  ^-bits  of  digital  information  are  represented  by  each 
signal.  The  possible  signals  to  be  transmitted  are  of  the  form 

Sm{t)  =  cos  2Trfmt  (2.1) 


where  E,  is  the  signal  energy  (E^  =  Jq  sl^{t)dt)  and  /„  is  one  of  the  M  frequencies 
used  to  distinguish  among  signals  where 


f  -  1  4. 

J  ^  rjy  *  ' 


Here.  Aj  is  an  integer  representing  the  spacing  (relative  to  1  /T)  between  each  of  the 
M  signals.  This  signal  is  modulated  with  an  appropriate  carrier  frequency  which 
is  removed  at  the  receiver. 

The  receiver  used  in  this  analysis  is  depicted  in  Figure  2.1.  After  removing  the 
carrier  frequency  and  low-pass  filtering  (with  cutoff  frequency  W  >  fmax)^  the  signal 
is  then  sampled  at  a  rate  N/T  {N  is  the  size  of  the  FFT)  so  that  the  outputs  of  the 
FFT  correspond  to  a  frequency  spacing  of  l/T.  With  this  sampling  rate,  only  the 
bins  with  a  spacing  Aj  from  Equation  2.2  will  be  used  in  the  decision  stage. 

The  input  to  the  receiver  is  the  transmitted  signal  corrupted  by  zero-mean, 
Gaussian-distributed  white  noise  T]{t)  with  power  spectral  density  (p.s.d.)  of  No/2. 
The  signal  may  also  experience  a  frequency  offset  of  /'  caused  by  doppler  effects, 
receiver  oscillator  drift  or  other  sources.  In  addition,  the  signal  experiences  slow 
Rician  fading  causing  the  amplitude  o  to  be  a  random  variable  following  a  Rician 


Figure  2.1:  With  a  FFT  receiver,  the  incoming  signal  is  sampled,  a  window  is 
applied,  and  the  magnitude  of  the  FFT  output  is  used  to  determine  which  signal  is 
received. 


distribution.  Rician-distributed  random  variables  in  signal  analysis  are  characterized 
by  the  ratio  of  the  power  of  the  direct  signal  component  to  the  power  of  the 
diffused  (or  faded)  component  2cr| 


/a(a)  = 


^exp 


2<t} 


where  /o(i)  is  the  modified  Bessel  function 


(2.3) 


Io(T)  ^  E 


.2n 


n=0 


22"(n!)2' 


(2.4) 


At  the  receiver's  input. 


r(t)  =  >/^acos[27r(/c  +  /„,  +  f')t]  +  Tj(t) 


(2.5) 


and  after  removing  the  carrier 


f(t)  =  \/2E,a  cos  [2;r(/„,  +  f')t]  +  ■q{t). 


(2.6) 


To  simplify  the  signal  representation,  the  following  substitutions  are  made; 

(2.7) 

(2.8) 


q  =  1  +  A/(?7T  —  1) 

£  =  fr 
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so  that  Equation  2.6  becomes 


r{t)  =  y/^sacos  +n{t)  (2.9) 

where  e  now  denotes  the  fractional  (of  l/T)  frequency  offset  and  q  represents  the 
multiple  of  l/T  of  the  sent  signal.  With  a  sampling  rate  of  N/T,  the  discrete  form 
of  Equation  2.9  is  represented  by 


r(n)  =  y^2£’,acos 


(2.10) 


As  seen  in  the  receiver  of  Figure  2.1,  the  discrete  signal  is  windowed  using  a 
discrete  windowing  function  w(n)  and  the  EFT  is  taken.  The  outputs  of  the  EFT 
can  be  expressed  in  real  (denoted  XR{k))  and  imaginary  A/(A’)  parts 


A'/?(A)  =  rin)w{n)cos 


2Ttkn 


V'' /  /  \  27r(^  +  e)n  27rkn  ,  .  ,  .  27rAn  1 

2^  <  w{n)J2E,a  cos - — - cos  — — - f-  7?(n)u»(n)  cos  — ^  >  (2.11) 

n=0  I  J\  N  A  I 


A'’/(A:)  =  ^  r(n)u’(n)sin 


2irkn 


V'  /  ^  /or"  27r(9  +  e)n  .  2Tkn  ,  ,  /  ,  .  27rA??  ] 

~  ^  - iv - - vi^)w{n)sm ^  .  (2.12) 


Keeping  in  mind  that  the  channel  experiences  slow  Rician  fading  so  that  a  is 
assumed  to  be  a  constant  over  the  period  of  interest  0  <  f  <  T,  and,  because  the 
received  signal  is  being  corrupted  by  random  noise  following  a  Gaussian  distribution. 
Ah(^)  3.r>d  Xj{k)  are  Gaussian  distributed  as  w'ell,  since  they  are  derived  from  linear 
combinations  of  Gaussian  processes.  As  a  result,  the  output  of  the  EFT  generator  is 
a  complex  Gaussian  process  of  the  form 


X(k)  =  XR{k)  +  jX,{k). 


(2.i;i) 
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The  necessary  variable  for  the  decision  criteria  will  be  the  magnitude  of  the 
outputs  of  the  FFT 

dk  ^  \Xik)\  =  ^/xUk)  +  X]ik).  (2.14) 

This  is  the  envelope  of  the  the  complex  signal  X{k).  As  stated  in  [3],  for  a  random 
complex  signal  X{k),  the  distribution  of  the  envelope  is 

1 


f{dk)  =  -^exp 


r  (  dkll 

^0  I 


where  (overbar  notation  denotes  expected  value), 


^  =  \X{k)\ 


r)l  = 


+  Xi(k) 


(2.15) 


(2.16) 


and 


^.V  = 


A-(A-)-A'(A:)|  J.  (2.17) 

The  mean  of  A(A')  (from  Equations  2.11  and  2.12)  for  a  fixed  amplitude  a  are 
easily  found  by  observing  that  /?(n)  is  the  only  random  quantity  and  is  zero-mean, 
causing  the  second  term  of  the  sum  to  go  to  zero  leaving  the  first  term  of  the  sum, 


N-l 


Xr  =  w(n)  cos 


2‘xn{q-\-c)  27rnk 
cos 


n=0 


N 


X 


A 

=  rrtR 


(2.18) 


^  ^  •  27rn(g  +  e)  .  27rnk 

A I  =  ay  2Ls  ^  u’(n}  Sin - — - sin 


n=0 


A 

=  m/. 


(2.19) 


The  variance  in  Equation  2.17  then  becomes 


^  2irnk  ,  ,  ,  ,  .  27rnA- 

Y  w{n)7]{n)cos--~  +  Y  ■w{n)T]{n)sm 

n=0  n=0 


A’ 
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(2.20) 


which,  because  of  white  noise,  the  autocorrelation  function  of  T){n)  is  given  by 

,No 


Rr,{T},p)  =  2W  -  p) 


(2.21) 
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where,  again,  W  is  the  cut-ofF  frequency  of  the  low-pass  filter  at  the  receiver  and  No/2 


is  the  p.s.d.  of  the  additive  white  Gaussian  noise.  Equation  2.20  is  reduced  to 


(N-l 

n=0 


.  o27rnk  ,  .  ^27rnk 

~ir  +  ^  w^n)sm^—— 

•'''  n=0 


r/  AT  N-l 


(2.22) 


Keep  in  mind  that  Equations  2.18  and  2.19  are  still  functions  for  a  given  k,  m 
and  e.  By  factoring  out  a  from  mn  and  m/  given  in  Equations  2.18  and  2.19,  and 
substituting  into  Equation  2.16, 


-f  m] 


=  a^{m%  ->r  mf) 


= 


(2.23) 


With  slow  Rician  fading.  Equation  2.15  is  actually  a  conditional  density  on  a  and 
becomes 

fd,[aidk\a)  =  ^exp  -~(4+a^/3k)  h  •  (2.24) 

The  density  function  of  dk  can  then  be  found  by  integrating  over  the  Rician- 
distributed  random  variable  a 

roc 

fdMk)  =  fdMk\a)fa{a)da 

=  r  f  (°4]  da 

Jo  cri-  k  a}  J  crj  [  cr?  I 


dk  ^  ( dl 

J-T^xp  -o  —  +  "T 
2  \(7l  a) 


r  ■“  r  adky/Wk\  ,  aA\  , 

Jo  \<^x  \^) 


(2.25) 


By  making  the  substitutions  Iq{x)  =  Jq{jx), 

a  ^-v 

2a\  2a]  2a\a] 


(2.26) 
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and  using  the  integral  (from  [4]) 


f  J^{ax)J^ibx)dx  =  ^r^exp  - 

Jo  2Q^ 


f  ab 


(2.27) 


gives 


f^Ad>.)  = 


\2Q^ 


1  \  + 
77^  - 0^2  ^2 


2ajc(rj 


di^k  .  21^^  1  1  ,  (-dkAsflTk 


(2.28) 


From  Equation  2.4,  it  is  seen  that  Iq{x)  is  an  even  function.  Replacing  the  as 
defined  in  Equation  2.26  into  Equation  2.28  results  in 

f  (A  \  ^  A'^a\ 

*  \0k(y'j  +  (7\  )  2<Tj^aj{l3kCrj  +  cr\) _ 

^  {  dkAs/^k  \ 

^  ^0  -a  .2 -,  -77  (2-29) 


X  Iq 


dkAy/^k 

0k(T)  + 


which  simplifies  to 


r  lA^-  (  ^  1  {dl  +  8kA^ 

Note  that  by  making  the  substitutions 


dkA\/Wk  N 

0k<y)  +  <^.Y  / 


(2.30) 


<^1  =  8k<T]  +  al 

Ok  =  AyJ^k-, 


(2.31) 


Equation  2.30  is  seen  to  be  Rician-distributed  since 

r  ^  dk  (  \di-\-  al\  ^  ( dkOLk 


(2.32) 


A  similar  result  was  found  in  [1]  where  the  parameters  al  and  q*.  are  now  altered  as 
a  result  of  Rician  fading.  Note  that  if  no  fading  is  present  (2cry  =  0  and  A  -  1),  then 
Equation  2.32  is  the  solution  found  in  [Ij. 
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III.  ERROR  ANALYSIS 


Now  having  an  expression  for  the  distribution  of  the  magnitudes  of  the  FFT 
output,  an  error  occurs  when  the  desired  frequency  component  (located  at  bin  q  of 
the  FFT)  is  less  than  any  of  the  other  FFT  bins  of  interest.  With  a  sampling  rate  of 
N/T,  the  bins  of  interest  are  those  corresponding  to  a  frequency  integer  multiple  of 
A/.  In  addition,  the  FFT  displays  symmetry  about  the  N/2  frequency  bin  so  only 
the  first  N/2  outputs  are  of  consequence.  Finally,  the  first  bin  corresponding  to  a 
d.c.  component  is  disregarded,  as  well  as  the  bin  located  at  N/2,  since  a  frequency 
component  in  this  bin  will  not  satisfy  the  Nyquist  criteria. 

As  further  discussed  in  [1]  (referencing  [5]),  the  union  bound  solution  is  carried 
out  for  the  probability  Ps  of  symbol  error.  From  Equations  2.18  and  2.19,  the  distri¬ 
bution  of  dk  will  differ  with  a  positive  or  negative  frequency  offset  e.  It  is  assumed 
that  a  positive  or  negative  frequency  offset  is  equally  likely.  In  addition,  each  of  the 
M  signals  are  assumed  equally  likely  to  occur  with  probability  l/M.  An  error  occurs 
when  the  magnitude  dk  of  a  FFT  bin  is  greater  than  the  magnitude  d,  (refering  to 
Equation  2.7)  of  the  bin  of  interest.  The  probability  of  symbol  error  is  then  given  by 


H  {Pr[d{q)  <  d{k)\  + c,q]  + Pr[d{q)  <  d{k)\  -  c.q]] 

ffc.  ^€1  +  A^(m  —  1) 

<  m  =1, 2,  .A/ 

(3.1) 

where,  from  [5].  for  Rician-distributed  random  variables 


Pr[d{q)  <  dik)]  = 

Q{y/d,y/b)  -  exp  /o(v^) 


(3.2) 
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with  the  following  definitions  made 


A 

a  = 


cm 


a: 


c  i  ^ 

Evaluation  of  the  Marcum-Q  function  is  described  in  Appendix  A.  Also,  the  depen¬ 
dency  of  e,  q  and  k  is  found  in  the  /3fc  term  in  the  definitions  made  in  Equation  2.23 
(with  reference  to  Equations  2.18  and  2.19). 

Parameters  for  analysis  include  the  direct-to-fading  ratio 


DTF  = 


2^2 


(3.3) 


and  the  signal-to-noise  ratio  at  the  receiver 


5,Vfl  S 

^o  No 

By  solving  each  of  the  above  equations  for  A^.  it  is  found  that 


(3.4) 


SNR 


A 

=  7 


(3.5) 


No  2{DTF  +  l) 

which  is  the  diffused  signal-to-noise  density  ratio.  Using  this,  making  the  substitutions 
from  Equation  2.31,  and  normalizing  the  transmitted  energy  (E,  =  1  in  Equation  2.18 
and  2.19),  then  a,  b  and  c  can  be  expressed  in  terms  of  the  given  parameters 

2DTF0k 


a  ~ 


(3.6) 


b  = 


2DTFI3. 


+  +  u-2(70 


(3.) 


c  = 


20,1  +  w  w^n) 


(3.8) 
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A  bound  on  bit  error  is  then  found  by 


Bit  Error  Rate 


IV.  WINDOWING  RESULTS 


Graphs  of  the  probability  of  bit  error  as  a  function  of  the  received  signal-to- 
noise  ratio  (per  bit)  were  constructed  for  various  direct-to- fading  ratios  .4/2cr^  for 
the  rectangular  and  Hamming  windowing  cases.  These  results  were  generated  from 
Equations  3.1  and  3.2  with  evaluation  of  the  Marcum-Q  function  and  the  modified 
Bessel  function  described  in  Appendix  A.  A  plot  with  low  fading  channel  (A^/2cr^  = 
100)  is  given  in  Figure  4.1  below.  The  case  studied  in  [1]  (A/  =  2.  M  =  8)  is  used  so 
that  a  direct  comparison  of  results  can  be  made.  In  the  ca.se  of  low  fading,  the  results 
are  comparable. 


(a)  Rectangular  (a)  Hamming 

Figure  4.1:  Probability  of  bit  error  for  low  fading  channel  (a)  rectangular,  and  (b) 
Hamming  windowing.  32-point  FFT.  M  =  8,  A/  =  2. 
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From  Figure  4.1  the  advantages  of  the  non-uniform  windowing  can  be  seen  as 
the  frequency  offset  becomes  larger  for  a  fixed  /„  and  T.  In  both  cases,  as  e  grows 
linearly,  the  energy  necessary  to  maintain  a  given  bit  error  rate  grows  exponentially. 
However,  the  advantage  of  the  Hamming  window  is  reflected  in  a  smaller  exponential 
increase  than  the  rectangular  window.  With  lower  values  of  frequency  offset  the 
signal  energy  necessary  is  larger  than  with  uniform  windowing,  but,  as  t  increases, 
the  advantages  of  non-uniform  windowing  become  apparent. 

A  visual  represenation  of  this  w'as  developed  in  [1]  by  graphing  the  amount  of 
SNR  degradation  necessary  to  maintain  a  fixed  probability  of  bit  error  as  a  function 
of  increasing  e  and  is  repeated  here.  Figure  4.2  gives  the  result  for  the  low  fading 
case  wdth  a  fixed  HER  of  10“^  (from  Figure  4.1). 


Rectangular 


0.5  ■  "  ' 

^  I^iycyo-Fading  Ratio  =  100.0  j 

%  at  ’  ’~0.2  OJ  O.F  0.5 

Fractional  Frequency  Offset  (t) 


Figure  4.2:  Trade-off  of  uniform  and  non-uniform  windowing  for  bit  error  rate  of 
10“^.  Non-uniform  windowing  becomes  more  advantageous  as  t  increases. 
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Similar  graphs  of  bit  error  rate  as  a  function  of  SNR  are  given  for  the  case 
of  a  slightly  diffused  channel  [A^ I2a^  =  3.16)  in  Figure  4.3.  The  results  display¬ 
ing  SNR  degradation  are  given  as  well.  From  Figure  4.3  we  see  that  with  a  more 
diffused  signal,  the  cross-over  fractional  frequency  offset  decreases.  This  phenomena 
was  studied  closer  by  generating  additional  results  in  the  same  manner  for  various 
direct-to-diffused  ratios  and  plotting  the  frequency  offset  cross-over  point  of  each  in 
Figure  4.4.  From  this  graph  we  see  that  the  jusfication  to  use  non-uniform  window¬ 
ing  is  very  sensitive  in  the  area  of  equal  direct  and  diffused  components.  However, 
as  the  direct  component  approaches  zero  (Rayleigh  fading),  non-uniform  windowing 
becomes  more  advantageous  when  minimal  frequency  offset  is  expected. 
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Figure  4.3:  Probability  of  bit  error  for  fading  channel.  Af2aj  =  3.16  with  (a) 
rectangular  and  (b)  Hamming  windowing.  SNR  degradation  for  a  fixed  BER  of  lO”" 
given  in  (c). 
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Figure  4.4:  Trade-ofF  for  various  direct-to-fading  ratios  ,4/2(7^.  As  the  channel 
approaches  Rayleigh  fading,  the  frequency  offset  necessary  to  justify  non-uniform 
windowing  decreases  (maintaining  a  BER  of  10~^). 


V.  ARQ  PROTOCOLS  AND  THEIR 
PERFORMANCE 

As  seen  in  the  previous  section,  noise  on  a  communication  channel  can  cause 
errors  in  the  reception  of  signals.  When  an  error  occurs  in  computer  communications, 
the  receiving  DLC  must  issue  a  command  to  the  sender  for  retransmission  of  the  frame. 
This  negative  acknowledgement  (NAK)  of  the  received  frame,  or  acknowledgement 
(ACK)  if  no  errors  occured,  can  be  a  relatively  short  frame  with  a  minimal  chance 
of  error  or  imbedded  in  a  longer  information  frame  with  a  higher  probability  of  error 
occuring  in  the  frame.  As  a  safeguard,  the  original  sending  node  assigns  a  maximum 
allowable  waiting  time,  or  timeout,  where,  upon  its  expiration,  the  frame  is  assumed 
to  have  been  received  in  error  and  is  then  retransmitted. 

The  efficiency  of  ARQ  schemes  is  dependent  upon  the  complexity  of  the  DLC  to 
implement  it;  the  more  complex  the  implementation,  the  higher  the  throughput  rate. 
Current  ARQ  protocols — stop-and-wait,  go-back-N  and  selective- repeat — are  based 
on  a  fixed  frame  length.  As  will  be  seen  from  the  throughput  curves  generated  for 
each  protocol,  a  higher  performance  may  be  achieved  in  areas  of  different  bit  error 
rates  with  different  frame  lengths. 

A.  STOP-AND-WAIT  ARQ 

The  stop-and-wait  ARQ  protocol  is  the  simplest  protocol  to  implement  and 
serves  as  a  good  example  to  illustrate  ARQ  schemes.  In  the  stop-and-wait  protocol, 
each  frame  is  transmitted  and  the  sending  DLC  then  stands  idle  waiting  for  a  response. 
If  a  negative  response  is  made  (either  by  a  NAK  or  by  a  timeout),  the  frame  is 
retransmitted.  This  process  continues  until  all  frames  have  been  transmitted.  As  an 


21 


example,  a  communication  is  illustrated  in  Figure  5.1  between  a  sending  DLC  (node 
A)  and  a  receiving  DLC  (node  B).  Here,  frame  ‘0’  is  sent  and  received  correctly  by 
node  B  who  ACKs  the  packet  by  requesting  the  next  frame  (frame  ‘1’).  An  error  was 
encounterd  with  frame  ‘T  and  node  B  NAKs  (negatively  acknowledges)  the  frame 
by  re- requesting  frame  ‘1’  which  is  retransmitted  by  node  A.  The  ACK  from  node  B 
encounters  an  error  and  at  the  conclusion  of  the  timeout,  frame  A  retransmits  the 
frame. 


*tranj 


Figure  5.1:  Stop-and-wait  protocol.  Frame  ‘0’  is  sent  and  ACKed  by  node  B.  An 
error  occurs  with  frame  T’  and  is  NAKed  by  B.  The  ACK  from  node  B  encounters 
an  error,  and  at  the  conclusion  of  the  timeout,  frame  A  retransmits  the  frame. 


1.  Performance  Analysis  Background 

Performance  for  each  of  the  three  protocols  will  be  developed  in  their  re¬ 
spective  sections.  Before  proceeding  with  the  analysis  for  the  stop-and-wait  protocol, 
the  following  background  used  in  all  three  performance  analyses  is  developed. 


a.  Assumptions 

The  following  assumptions  are  made  throughout  the  performance  analysis: 

•  Node  (A)  continuously  transmits  to  node  (B).  —  With  node  A  continuously 
transmitting,  there  are  no  queuing  delays.  Because  node  B  has  no  information 
to  send  to  node  A,  acknowledgement  frames  will  only  contain  the  ACK  or  NAK. 

•  An  acknowledgement  is  always  received  and  is  never  received  in  error.  —  This 
assumption  can  be  made  because  of  the  relatively  short  length  of  the  acknowd- 
edgement  frame. 

•  The  propagation  delay  is  known.  —  This  stipulation  is  only  necessary  to  es¬ 
tablish  an  appropriate  window  size  with  the  go-back- and  selective  repeat 
protocols. 

•  Processing  time  at  each  node  is  negligible  and  can  be  considered  included  in 
the  propagation  delay. 

•  Bit  error  is  independently  and  identically  distributed  for  all  bits  within  a  frame. 

b.  Average  Transmission  Attempts 

If  the  probability  f*;.  o^  bit  error  is  independent  and  identically  distributed  for 
each  bit  within  a  frame  {L  bits  long),  the  probability  of  a  frame  being  received  in 
error  is  given  by 

Pe  =  \-  Psuccss  =  1  -  (1  -  (5.1  ) 

The  probability  that  it  will  take  i  attempts  to  successfully  transmit  a  packet  can  be 
expressed  as 

Pr[number  of  attempts  =  i]  =  —  Pe).  (5.2) 


2.3 


The  expected  number  of  attempts  can  be  computed  by  ’ 


c.  Throughput  Definition 

Throughput  is  defined  as 

A  total  bits  successfully  transmitted  in  a  given  time  interval  ('6) 

^  bit  capacity  of  the  channel  in  the  same  interval 

This  can  also  be  expressed  as  the  ratio  of  the  time  to  transmit  a  frame  t^rans  to  the 

average  time  to  transmit  a  frame  without  error  t. 

Because  the  information  packet  is  appended  with  overhead  bits  L  some  of  which 
contains  coding  information,  the  throughput  is  scaled  by  this  ‘coding'  rate.  This 
results  in  a  throughput  efficiency  given  by 

A  itrans  (T  () 

P  =  - - (o.O 

2.  Stop-and-Wait  Performance 

The  total  time  the  transmitter  expends  for  each  frame  is  a  sum  of  the  trans¬ 
mission  time  of  the  frame,  transmission  time  of  the  acknowledgement  tack  and  two 
propagation  delays  tprop  (which  are  assumed  to  contain  the  processing  delay)  for  the 
round-trip  cycle.  This  gives  a  total  time  for  each  transimission  attempt  of 


(1  -  p,)E 

oo 

ii-p,)j:zpr^ 

1 

i-p. 


t=i 


(5.3) 

(5.4) 

(5.5) 


Pl  —  ^trans  d"  T  tack- 


(5.8) 


'The  infinite  sum  is  shown  to  be  convergent  by  taking  the  derivative  of  the  geometric  infinite 
series  solution  where  x  <  1 


1 


dx 


n  =  0 


dx  l  —  X 
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Using  the  results  of  Equation  5.5  for  the  average  number  of  attempts,  the  average 
time  to  transmit  a  frame  succussfully  is  given  by 


taw  — 


"b  "b  tack 


1-Pe 


(5.9) 


By  applying  the  definition  of  throughput  given  in  Equation  5.7  and  defining  the 
number  of  frames  that  can  be  sent  in  one  round-trip  propagation  delay  cis 


A  2i 
a  =  — 


prop 


"b  tuck 


t 


trans 


(5.10) 


the  throughput  of  the  stop-and-wait  strategy  is 


Paw  — 


(5.11) 


l-|-a  L  l-j-a  L 

Theoretical  curves  for  Equation  5.11  are  given  in  Figure  5.2  for  various  frame  lengths. 


The  simulations  described  in  Appendix  B  were  performed  and  the  results  are  given 
here  as  well. 

Although  the  stop-and-wait  protocol  is  easy  to  implement,  the  efficiency  of  the 
system  is  compromised  by  waiting  for  a  response  after  each  packet  is  transmitted.  The 
following  two  protocols  do  not  stand  idle  waiting  for  a  response  and,  subsequently, 
have  better  expected  throughput  than  the  stop-and-wait. 


B.  GO-BACK-A^ 

In  the  go-back-jV  protocol,  node  A  continuously  sends  frames  to  node  B.  Node 
B  keeps  node  A  aware  of  its  status  by  responding  with  an  acknowledgement  which 
contains  an  RN  (received  number)  indicating  that  all  frames  preceding  RN  have  been 
received  correctly.  RN  also  indicates  that  all  frames  from  RN  onward  need  to  be  sent 
by  A.  As  frames  are  received  without  error,  RN  is  incremented. 

In  the  event  of  an  error,  RN  will  remain  constant  as  illustrated  in  Figure  5.3. 
There  is  a  limit  to  the  number  of  frames  node  A  can  send  beyond  RN  so  that  a  frame 
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in  error  can  be  detected  and  retransmitted  by  node  A.  The  number  of  frames  that 
can  be  transmitted  beyond  RN  is  the  A’ in  go-back- A'  and  is  referred  to  as  the  window- 
length.  Because  RN  remains  constant^ when  an  error  occurs,  node  A  will  stop  as  soon 
as  the  end  of  the  window  is  reached  (5 A’  -f  1  =  RS  +  A’)  and  go  back  A'  frames 
and  begin  retransmission  of  frame  RN  and  all  subsequent  frames.  To  minimize  the 
number  of  frames  needed  to  be  retransmitted  in  the  event  of  an  error,  the  value  of  A’ 
is  set  to  the  number  of  frames  that  can  be  transmitted  in  the  time  interval  given  in 
Equation  5.8.  Using  Equation  5.10, 


A'  =  a  -I-  1 . 


(5.12) 


With  the  go-back-A’  protocol,  the  transmitter  can  send  frames  continuouslv. 
However,  upon  each  frame  received  in  error,  the  protocol  requires  that  the  .sending 
DLC  go  back  A’  frames  and  re.send  them.  Thus,  when  successful  on  the  /th  attempt. 
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Go- back-4 


i 


node  B 

Figure  5.3:  Go-back-A  protocol  for  N  =  4.  With  frame  2  received  in  error.  RN 
stays  constant  at  2  until  node  A  has  reached  the  end  of  the  window.  Node  .N  detects 
the  error  and  must  ”go  back  4’'  frames. 

the  DLC  has  sent  a  total  of  1  +  {?  —  1)A^  frames.  The  average  number  of  attempts 
was  given  in  Equation  5.5  which  results  in  an  average  number  of  frames  (/)  that  need 
to  be  sent  of 


/  = 

l+Pe(JV-l) 

l-Pr  ' 


(5.13) 


Using  Equation  5.12,  the  average  time  to  transmit  a  frame  successfullv  is 


^  —  ^transf 

Uransi  I  +  G )  ,  -  ,  , 

1  -  I\  ■ 

Applying  the  definition  given  by  Equation  5.7,  the  theoretical  througliput  for  the 
go-back- Af  protocol  becomes 


1  -f  aP,  L 


(5.15) 
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where  the  dependency  of  Pc  on  L  is  shown  in  Equation  5.1.  A  graph  of  theoretical 
throughputs  for  various  frame  lengths  is  given  in  Figure  5.4  with  the  simulation 
results  superimposed.  It  can  be  seen  from  the  graph  that  an  adaptive  strategy  will 
not  improve  the  performance  of  the  go-back- protocol.  A  frame  length  of  1024  could 
be  used  at  all  bit  error  rates  and  only  a  fractional  percent  of  throughput  would  be 
lost  at  higher  bit  error  rates. 
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Figure  5.4:  Performance  of  go-back- A'  for  various  frame  lengths.  A  frame  length 
of  1024  could  be  used  at  all  bit  error  rates  with  only  minimal  loss  of  performance  at 
higher  BER. 


C.  SELECTIVE-REPEAT  ARQ 

The  selective- repeat  protocol  is  similar  to  the  go-back- A^  e.xcept  that  when  an 
error  occurs,  only  the  frame  that  is  in  error  needs  to  be  retransmitted.  Here  it 
is  important  to  acknowledge  each  frame  (both  positive  and  negative)  so  that  the 
transmitter  accounts  for  all  frames  sent  out.  If  a  response  is  not  received,  that  frame 
reaches  its  timeout  and  is  retransmitted. 


Frame  5  timeout 


Figure  5.5;  Selective  repeat  protocol. 


The  obvious  advantage  of  selective- repeat  over  go-back- A"  is  that  only  the  frame 
in  error  is  retransmitted  as  opposed  to  sending  all  following  N  —  1  frames  in  the 
go-back- A^  protocol.  However,  this  requires  that  the  receiving  DLC  have  adequate 
memory  to  store  frames  that  it  receives  and  to  be  able  to  sort  them  as  they  become 
available.  If  memory  in  the  receiver  is  full,  then  all  incoming  frames  are  NAKed. 
regardless  of  any  errors. 

The  throughput  analysis  for  selective-repeat  (with  an  infinite  receiver  buffer 
size)  is  straightforward  from  Equation  5.5.  Because  the  only  delay  experienced  is 
that  of  the  retransmission  ol  a  single  frame  when  an  error  occurs. 


t  = 


f  irons 

l-Pe 


(5.16) 


which  results  in  a  theoretical  throughput  for  selective- repeat  with  an  infinite  receiver 
buffer  as 

p„  =  {l- (.5.17) 

Plots  of  theoretical  throughputs  with  infinite  receiver  buffer  length  are  given  in  Fig¬ 
ure  5.6  with  the  results  obtained  in  the  simulation. 


Q  =  I- 

3  =  \- 

1  =  \-pt- 


30 


f- 


Ir- 

t- 
0.9  r 

I 

I 

0.8 1- 

1 

0.7  H 


3  0.6  h 

^  I 

on 

3 


0.5 


) 

0.4- 

I 

0..7‘- 

I 

0.2  r 

I 

0.11- 


0- 

10' 


Lower  Bound  (22K  Buffer) 


Simulation  (64K  Buffer) 


■■  /'  i  ■  Infinite  Buffer 


C= 32000  b/s 
prop =700  msec  , 

ACK=32  bits  J 

over=32  bits 


10’  W 

Probability  of  Bit  Error 


10^ 


10  2 


Figure  5.7:  Performance  of  SR  with  finite  receiver  buffer  (fixed  length  L  =  1024). 
N  =  22  is  the  number  of  frames  transmitted  in  a  round-trip  delay.  From  the  analysis 
of  [6],  a  lower  bound  on  throughput  is  given. 


A  comparison  of  the  upper  bound  (infinite  receiver  buffer)  and  the  lower  bound  given 
in  Equation  5.18  is  given  in  Figure  5.7  for  a  fixed  length  of  I  =  1024.  For  the 
simulation,  a  receiver  buffer  size  of  64K  is  used  instead  of  restricting  it  to  .V  (.V  =  22 
for  the  parameters  of  the  simulation).  A  buffer  size  of  64K  was  believed  not  to  be 
too  excessive,  yet  imposes  some  restriction  on  the  system. 
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VI.  ADAPTIVE  ARQ  PROTOCOLS 


Implementation  of  the  previously  described  ARQ  schemes  remains  the  same  in 
an  adaptive  environment.  The  adaptation  takes  the  form  of  varying  the  length  L 
of  the  frames.  With  an  increased  bit  error  rate,  the  increase  in  frame  error  rate  in 
Equation  5.1  can  be  offset  by  decreasing  the  length  of  the  frame.  Alternatively,  if 
the  bit  error  rate  of  the  channel  decreases,  a  longer  frame  length  will  decrease  the 
percentage  of  a  frame  dedicated  to  overhead  (assuming  the  overhead  remains  constant 
for  varying  lengths). 


A.  ADAPTIVE  STOP-AND-WAIT  USING  A  SEMAPHORE 


The  decision  to  change  the  length  of  the  frame  was  made  in  three  different 
schemes.  The  first  scheme  is  taken  from  [7]  and  invokes  a  counting  semaphore  m 
that  is  incremented  each  time  a  frame  is  acknowledged  and  decreased  each  time  a 
frame  needs  to  be  retransmitted.  The  length  of  the  frame  is  then  dependent  upon 
the  current  value  of  the  semaphore  by 


'I  7o  <  m  <  7i 

2/  7i  <  m  <  72 


I  nl  7„_i  <Tn<^n- 


(6.1) 


For  this  thesis,  I  =  128  and  n  =  4. 

Figure  6.1  shows  the  gain  associated  with  an  adaptive  stop-and-wait  protocol. 
This  scheme  adds  little  complexity  to  the  DLC  implementation  and  provides  the 
desired  results. 
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L=1024 


L=512 


L=256 


L=128 


C=  32000  b/s 
prop  =  20  msec 

ACK=32  bits 
over  =  32  bits 

\  —  Fixed  L 

\  -o-  Adaptive 


\,  Adaptive 


(Simulation  Results) 


Probability  of  Bit  Error 

Figure  6.1;  Adaptive  stop-and- wait  performance.  Improved  performance  is  obtained 
in  the  region  of  high  BER. 

B.  ADAPTIVE  SELECTIVE-REPEAT  USING  THRESHOLDS 

.Another  criteria  compared  the  number  r  of  frames  transmitted  after  a  fixed 
\alue  b  of  frames  were  received  in  error  to  a  threshold  to  determine  if  the  frame  length 
should  be  increased  (doubled)  or  decreased  (halfed).  Two  schemes  were  implemented 


where  either  a  single  threshold 


or  a  dual  threshold 


L/2  6/r  >  7 

2L  b/r  <  7 


^/2  Pl/i  >  pL 

L  =  <  2L  P2l>  PL  (6.3) 

L  otherwise 

was  used.  The  computation  of  pi/2'  PL  and  P2L  are  based  on  the  experienced  bit  error 
rate.  The  experienced  P\,  can  be  computed  from  the  ratio  of  transmitted  frames  to 
frames  received  in  error  by  applying  the  inverse  of  Equation  h.\ 


=  1  -  P,  =(1  -  l\f 
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1/Z, 


P6  =  l- 


(6.4) 


Throughputs  for  lengths  2L  and  T/2  can  then  be  computed  for  the  selective  repeat 
protocol  with  Equation  5.17.  The  decision  criteria  of  Equation  6.3  then  simplifies  to 


L  =  i 


1/2  1  -e/iL-£)>[ir~b)/r]^^^ 
2L  1  -  e/{2L  -£)<{r-  b)/r 
L  otherwise 


(6.5) 


Table  6.1  can  be  constructed  for  the  dual  threshold  selective-repeat  case.  The  deter¬ 
mination  of  7  in  Equation  6.2  was  taken  as  an  average  (7  =  0.88)  of  the  thresholds 
necessarv  to  cause  transition  from  table  6.1. 


L-^2L 

L  1/2 

L  =  128 

0.857 

N/A 

L  =  256 

0.933 

0.7.35 

L  =  512 

0.968 

0.871 

L  =  1024 

N/A 

0.936 

TABLE  6.1:  Adaptive  SR  transition  levels.  The  necessary  value  of  (r  —  b)jr  to 
cause  transition  either  up  or  down  with  the  dual  threshold  implementation  [k  =  32 
bits). 

The  performance  of  the  adaptive  selective-repeat  for  a  fixed  buffer  length  is 
given  in  Figure  6.2.  Implementation  of  this  scheme  further  increases  the  complexity 
of  the  selective- repeat  protocol.  Now  the  transmitting  node  must  keep  track  of  the 
length  of  each  frame  it  is  sending  out.  In  addition,  re-arrangement  of  the  frames  at 
the  receiving  node  becomes  quite  complex.  For  example,  if  a  frame  is  N.'\Ked  and 
the  decision  to  change  frame  lengths  is  made,  the  packet  sent  in  error  will  now  be 
divided  between  two  frames  each  of  which  may  or  may  not  be  received  correctly  in 
order  upon  retransmission. 
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Probability  of  Bit  Error 


Figure  6.2:  Adaptive  selective- repeat.  Both  dual  and  single  threshold  results  are 
given.  Performance  of  the  single  threshold  implementation  is  slightly  lower  in  the 
areas  where  the  decision  to  change  lengths  can  ’’toggle”  between  two  lengths. 

C.  GO-BACK-AARQ 

The  go-back- iV  protocol  does  not  lend  itself  to  an  adaptive  environment  as  was 
seen  in  Figure  5.3.  The  difference  in  throughput  rates  between  the  largest  frame 
length  used  (1024)  and  the  maximum  achievable  throughput  is  negligible.  This  is 
because  the  inefficiency  of  go-back-rV  arises  from  the  necessity  to  go  back  N  frames 
when  an  error  in  a  frame  occurs.  This  is  depicted  in  Figure  6.3  where,  if  the  first 
frame  is  in  error,  the  X  —  1  following  frames  composed  of  L  bits  each  result  in  a 
possibly  unnecessary  retransmission  of  {N  —  \)L  bits.  If  the  frame  length  is  halved 
(and  X  doubled  to  satisfy  Equation  5.12),  the  probability  of  having  to  re-transmit 
the  same  number  of  bits  is  the  same  because  the  probability  of  transmitting  two 
successive  frames  of  length  Lj'I  correctly  is  equal  to  the  probability  of  transmitting 
a  single  frame  of  length  L  without  error.  Although  some  advantage  may  be  gained 


35 


in  successfully  transmitting  the  first  frame  of  length  L/2,  there  is  also  an  increase 
in  percentage  overhead  for  each  frame.  For  example,  with  a  BER  of  10''*,  an  initial 
length  of  I  =  1024  and  overhead  (  =  32,  if  the  frame  length  is  halved,  P,ucc  increases 
by  only  4.7%  while  the  percentage  of  frame  overhead  increases  3.1%.  This  results  in 
negligible  improvement  in  the  areas  of  higher  bit  error  rates. 


Repeated  if  proceeding  frame  in  error 


(1  -  =  P,.cc  =  (1  - 


Figure  6.3:  Inefficiency  of  the  go-back- .V  protocol  occurs  because  frames  following  an 
error  must  be  repeated.  If  the  frame  length  is  halved,  the  probability  of  retransmitting 
the  same  number  of  bits  remains  the  same. 


VII.  CONCLUSIONS 


In  an  M-ary  FSK  scheme  implemented  using  sampled  data  at  the  receiver  and  a 
FFT  to  determine  which  signal  was  sent,  it  was  found  that  non-uniform  windowing  can 
be  used  to  increase  the  performance  of  the  system  in  the  presence  of  large  frequency 
offsets.  Implementing  non-uniform  windowing  causes  greater  separation  between  fre¬ 
quency  components  of  the  FFT  when  frequency  shifts  occur  between  transmission 
and  reception.  Non-uniform  windowing  also  causes  attentuation  of  the  signal  making 
detectability  more  difficult.  A  balance  between  improved  frequency  separation  and 
attenuation  is  achieved  as  the  frequency  offset  becomes  more  dramatic.  In  the  case 
of  a  fading  channel,  the  frequency  offset  necessary  to  justify  non-uniform  windowing 
becomes  smaller  as  the  channel  approaches  Rayleigh  fading. 

The  application  of  an  adaptive  ARQ  protocol  was  shown  to  improve  the  through¬ 
put  of  the  stop-and-wait  and  selective-repeat  ARQ  protocols.  With  an  adaptive  frame 
length,  throughput  can  be  increased  by  decreasing  the  probability  of  frame  error  in 
high  bit  error  rate  situations  through  a  reduction  in  the  frame  length.  With  low 
bit  error  rates,  a  larger  frame  length  reduces  the  percentage  overhead  in  the  frame 
thereby  increasing  the  throughput.  The  go-back- protocol  does  not  lend  itself  to 
adaption.  From  theoretical  analysis  and  simulation,  it  was  found  that  a  fixed  frame 
length  can  be  used  in  the  go-back-A  protocol  through  all  bit  error  rates  with  minimal 
loss  of  performance. 
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APPENDIX  A:  ERROR  COMPUTATION 


The  probability  given  in  Equation  3.2  was  computed  using  the  enclosed  C  pro¬ 
gram  code.  Values  for  the  Marcum-Q  function  and  the  Modified  Bessel  Function  were 
calculated  to  within  a  designated  tolerance  (10“®)  using  power  series  expansions  of 
the  functions  from  [4]  and  [5] 


Q{c^^0)=< 


'  1 

e-^V2 

i(l  +  /o(a2))e-<^^ 


exp 

exp 


0^+c.^ 


2 

0^+a^ 


if  d  =  0 
if  Q  =  0 
if  Q  =  d 
if  a  <  /? 

In{a0)  if  Q  >  ;? 


(.4.1) 


where 


r  1  if  X  =  0 

1  (x/2r  ESo  IfgS  otherwise 


(A.2) 
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1.  Bit  Error  Program  Code 


/***^*^i*j^********t***************************************** 


*  Pefad.c  * 

*  Calculates  the  probability  of  bit  error  for  a  Rician  * 

*  fading  channel  and  a  DFT  based  receiver  for  MFSK  * 

*  * 

*  by:  Chris  G.  Kmiecik  ♦ 

♦  date:  22  April,  1990  * 


#include  <stdio.h> 

#include  <math.h> 

tdefine  pi  3.141592654 
#define  sqrt2  1.414213562 
#define  epsilon  0.00000001 

#define  N  32 

#define  M  8 

#define  1  3 

#define  U  16 

#define  df  2 

double  sindouCM]; 

mainO 

{ 

double  mipO,  mrp(),  I(),  Q(),  powerO,  rectangleO,  hammingO,  banning () ; 

double  Perror,  lastPerror,  Enoise; 

double  a,  b,  c,  niunb,  numc,  denom; 

double  meanrq,  meaniq,  betaq; 

double  ine2Lnrk,  meauiik,  betak; 

long  factO: 

FILE  *output ; 
int  n,  k,  q,  einc; 

float  SHRdb,  SHR,  DTF,  DTFinc,  ratio,  e,  freqe; 
systein("nn  Pef adoutHAMH.mat") ; 


Enoise=0.0;  /♦  */ 

for  (n=0  ;  n  <=  I-l  ;  n++')  {  /*  Compute  the  values  ♦/ 

window [n+1] =h2unming(n) ;  /*  for  the  window  at  */ 

Enoise=Enoise+W*window[n+l]*window[n+l]/2.0;  /♦  each  ’n’  and  the  energy  */ 

>  /*  */ 
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lor  (DTFinc=2.0  ;  DTFinc  >=  -2.0;  DTFinc  =  DTFinc-0 . 25){/*  Specify  Direct  to  ♦/ 
DTF=pow(10.0,DTFinc) ;  /♦  fading  ratio  */ 


(e=0.0  :  e  <=  0.5  ;  e=e+0.1)  { 

/♦ 

Cycle  through  various 

*/ 

/♦ 

frequency  shifts 

♦/ 

Perror=3.0;  lastPerror=4.0: 

SIRdb=2; 

while ((Perror  >=  0.00001) 

/♦ 

Continue  computing  Pb 

♦/ 

kk(SIRdb  <=  99) 

/♦ 

until  SIR  becomes  large 

*/ 

kkdastPerror-Perror  >=  .  00001)  ){ 

/* 

or  error  becomes  small 

♦/ 

lastPerror=Perror ; 

/♦ 

or  error  does  not  change 

♦  / 

SHRdb++ : 

SHR=l^pou(10.0,SHRdb/10.0): 

/♦ 

Convert  Eb  to  Es 

♦  / 

ratio=SNR/(DTF+l) ; 

Perror=0 . 0 ; 

for  (q=l  :  q  <=  M*df-1  ;  q=q+df)  { 

/* 

Each  value  of  q  donates 

*/ 

me2mrq=mrp ( q , e , q) ; 

/♦ 

to  Pe.  First 

♦/ 

meaniq=mip(q,e,q) ; 

/♦ 

find  the  me2ai  of  the  real*/ 

betaq=meanrq*meanrq+maaniq*meaniq; 

/* 

3Uid  imaginary  paurts  then 

*/ 

numc=0 . 5*betaq*ratio+Enoise ; 

/♦ 

specify  all  peirauneters 

♦/ 

numb=(SNR-ratio)*betaq; 

/♦ 

changing  with  q  only. 

♦/ 

for  (einc  =  1;  einc  <=  2;  einc++)  { 

/♦ 

Include  both  plus  amd 

*/ 

lreqe=e*poHer(-l .0,einc) ; 

/* 

minus  Doppler  shift 

♦/ 

for  (k=l  ;  k  <=  (H/2)-l  ;  k=k+dl) 

if  (k  !*  q)  { 

meeuu:k=mrp(q,freqe,k) ; 

/* 

Calculate  values  in 

*/ 

meanik=mip(q,freqe,k) ; 

/♦ 

equations  2.17,  2.18, 

*/ 

betak=meanrk*meanrk+meanik*me2mik;  /*  2.22,  3.1.  ♦/ 

c=nujBc/(0 . 5*betak*ratio+Enoise) ; 
denom=0 . 5* (betaq+betak) *rat io+2 . 0*Enoise ; 
a= ( SIR-r at io ) ♦bet  ak/ denom ; 
b=nuinb/denoni ; 

Perror=Perror+Q(sqrt(a) ,sqrt(b) )  /♦  Equation  3.2  ♦/ 

-(c^exp(-0.6^(a+b) )^I(0,sqrt(a*b) ) )/(c+l) ; 

}  /*  if  (k  !=  q)  ♦/ 

}  /♦  for  (k=l  ;  k  <=  (H/2)-l  ;  k=k+df)  ♦/ 

}  /♦  lor  (einc  =  1;  einc  <=  2;  einc++)  ♦/ 

>  /♦  for  (q=l  ;  q  <=  M^df-1  ;  q=q+dl)  */ 

Perror=Perror/(2.0^H) ;  /♦  Equation  3.1  */ 

if  (Perror  <=  0.2)  { 

output =f open ("PeladoutHAHM. mat", "a") ; 
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Iprintf  (output , "*/.10 . 81  */,4.21  */,5.2f\n"  .Perror  ,SKRdb,e,DTF) ; 

lclose(output) ; 

> 

}•  /*  uhile  (Perror  >=  0.00001)  */ 

}  /*  lor  (e=0.0:  e  <=  0.5  ;  e=e+0.1)  •/ 

>  /*  lor  (DTFinc=2  ;  DTFinc  >=  -2;  DTFinc  =  DTFinc  +  O.02)  */ 

}  /*  close  main  ♦/ 
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♦  * 

♦  Function  mrp  ♦ 

*  * 
^^^^i^i^^^^^^^^^^^^mi>^*■**^^*%****^^^**■l^’¥***********i^*********/ 

double  mrpCq.e.k) 
int  k,q: 

■'loat  e; 

< 

int  n : 

double  ans; 

ans=0.0; 

for  (n=0;  n  <=  I-l;  n++)  { 

ans=eais+cos(2*pi*n*(q+e+k)/H) ♦window Cn+l]+cos(2*pi*n*(q+e-k)/H)*HindoH[n+l3 ; 

> 

if  ((e  ==  0.0)  kft  (k  !=  q))  ans  =  0.0; 

ans=aLns/sqrt2 ; 
return (ans) ; 

> 

♦  ♦ 

•  Function  mip  ♦ 

♦  * 

double  mipCq.e.k) 
int  k,q; 

float  e; 

{ 

int  n; 

double  ans; 

euis=0 . 0 ; 

for  (n=0;  n  <=  N-1;  n++)  < 

ans=ans+s in(2*pi*n* (q+e+k)/N)*window[n+l]-sin(2*pi^n*(q+e-k)/N) ‘window [n+1] ; 

> 


if  ((e  ==  0.0)  kft  (k  !=  q))  ans  =  0.0; 


ans=ans/sqrt2 ; 
return (ans) ; 

} 
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* 

* 

* 

4i4>*4>*«4 


Marcum-Q  Function 


double  Q(alpha,beta) 

double  alpha,  beta; 

{ 

double  ans,  e,  inc; 

int  n; 


if  (beta  ==  0.0) 
ans  =  1.0; 

else  if  (alpha  ==  0.0) 

ans  =  exp(-0.5*beta*beta) ; 

else  if  (alpha  ==  beta) 

ans  =  0 . 5*(1 .0+I(0,alpha*alpha)*exp(-l*alpha*alpha) ) ; 

else  if  (alpha  <=  beta)  ( 
ans  =  0.0; 

e  =  exp(-0.5*((alpha*alpha)+(beta*beta))) ; 

n=0;  inc  =  1.0; 

while  (inc  >=  epsilon)  •( 

inc=e*power( alpha/bet a, n)*I(n,alpha*beta) ; 
ans  =  ans  +  inc; 
n++ ; 

} 

} 

else  { 

ans  =  0.0; 

e  =  exp(-0.5*((alpha*alpha)  +  (beta>fbeta))) ; 

n=l;  inc  =  1.0; 

while  (inc  >=  epsilon)  { 

inc=e*power ( alpha/beta , n) *1 (n , alpha*beta) ; 

ans  =  ans  +  inc ; 

n++; 

> 

ans=l . 0-ans; 

> 


return (ans) ; 

} 
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Modified  Bessel  Function 


double  I(v,x) 

int  V ; 

double  X ; 


double  ans,  inc; 

int  n; 


if  (x  ==  0.0) 
ans  =  1 : 

else  { 
x=x/2.0: 

ans=poHer(x,v)/fact((long)  v); 

inc=ams ;  n=l: 
dhile  (inc  >=  epsilon)  { 
inc=inc*x*x/(n*(v+n) ) ; 
ans=ans+inc ; 
n++ ; 

> 

> 


return(ans) ; 

} 
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♦  * 


♦  Fixnction  rectangulcir  ♦ 

♦  * 
Jt^t*****t***********i************t*mm******************/ 

double  rectangle (n) 
int  n; 


retumd  .0)  ; 

> 


«  * 

♦  Function  Hamming  * 

*  * 
m*********’¥*****m*************»‘*******m***************/ 


double  hamming (n) 
int  n; 


{ 

return(0 . S4-0 .46*cos (2*pi*n/(N-l) ) ) ; 

} 
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*  * 

♦  Factorial  ♦ 

4>  * 

long  lact(x) 
long  x; 


long  i ,  ems ; 

ans=l ; 
il  (x  ==  0) 
ans  =  1; 
else 

for  (i  =  1;  i  <=  x;  i++)  ans  =  ans*i; 

return(ans) ; 

> 

*  * 

*  POWER  ♦ 

♦  * 

double  power (x,n) 
int  n; 

double  x; 

{ 

double  p; 
int  recip: 

if  (n  <  0)  •( 
recip=l ; 
n=abs(n) ; 

} 


for  (p=1.0;  n  >  0;  — n)  p=p*x; 

if  (recip  ==  l)  p=1.0/p; 

return(p) ; 

} 
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APPENDIX  B:  ARQ  SIMULATIONS 

The  ARQ  simulations  were  carried  out  using  a  C-program.  Although  one  node 
would  be  transmitting  to  another  in  the  application  of  an  ARQ  protocol,  the  simu¬ 
lation  only  involves  one  node.  A  doubly-linked  list  maintains  a  record  of  jobs  to  be 
performed  in  chronological  order  by  the  node.  The  node  can  either  transmit  (’SEND  ) 
or  receive  an  acknowledgement  (’RECV’)  for  a  given  frame. 

For  each  transmission,  the  time  of  completion  of  transmission  and  the  computed 
tinre  of  reception  (based  on  Equation  5.8)  are  added  to  the  list  of  jobs.  .4t  tnc 
completion  of  transmission,  the  node  then  transmits  the  next  frame  if  the  protocol 
allows  it. 

At  the  execution  of  ’RECV',  a  uniformily  distributed,  pseudorandom  number 
between  zero  and  one  is  compared  to  the  theoretical  frame  error  rate  (Equation  5.1). 
If  the  comparison  dictates,  a  NAK  is  recorded  by  the  node  and  the  appropriate 
retransmission  protocol  is  applied.  Particulars  for  implementing  retransmission  in 
each  protocol  are  given  below. 

1.  Selective-Repeat 

This  simulation  served  as  the  template  for  the  others  (see  Section  B.4)  in  that  go- 
back- A' and  stop-and-wait  are  ’’subsets"  of  this  program  .  The  adaptive  nature  of  the 
program  is  contained  in  a  section  added  after  receiving  each  frame’s  acknowledgement, 
(see  Figure  B.l). 
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2.  Go-Back- 


Before  sending  a  frame,  the  window  is  checked  to  see  if  it  is  full.  If  the  window 
is  full  SN  is  reset  to  the  beginning  of  the  window  (Figure  B.2). 

3.  Stop-and-Wait 

Here,  only  the  time  to  receive  each  packet  is  loaded  onto  the  linked  list.  When 
a  frame  is  received,  it  is  either  ACKed  or  NAKed  and  a  frame  is  sent  out  (see  Fig¬ 
ure  B.3).  If  the  frame  is  ACKed,  the  number  of  successful  bits  is  incremented. 

4.  Selective-Repeat  Program  Code 


*  sr.c  ♦ 

*  * 

*  Simulation  of  selective  repeat  protocol  with  data  generated  * 

*  for  BER  from  lE-6  to  lE-2.  * 

* 

*  Receiver  buffer  size  is  kept  track  of  by  the  sender  * 

*  * 

*  by:  Chris  G.  Kmiecik,  LT,  USCG  * 

*  date:  05  Februaxy,  1990  * 


#include  <stdio.h> 

#include  <malloc.h> 

#include  <math.h> 

#define  HALLOC(x)  ((x  ♦)malloc(sizeof (x) ) ) 
tdefine  max_FH  4096 

#define  max_buffer  65636.0 

typedef  struct  job  { 

float  time,  length; 

chcir  typo  [5]; 

int  FN ; 

struct  job  ♦previous,  ♦following; 

>  job.type; 

float  logBER; 

job.type  ♦top.job; 

FILE  ♦of,  ♦tf,  ♦rf,  ♦ttf; 

main( ) 

{ 

int  framestatus[max_FH] ,f ramelength[max_FN] ; 
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int 

int 

float 

float 

float 

float 

float 

job_type 

unsigned  long  int 
char 


queue ; 

buffer,  begin_windos,  SI,  i; 
completed,  m,  Ps; 

BER,  simulation_length; 

capacity,  propogation,  overhead,  ack_length; 
rand,  length; 

8uccessful_bits ,  throughput,  tr,  tt; 

*temp,  eprevious.job,  *job,  eneBjobO; 
seed; 

estreat () ; 


seed=1234: 

systemC'rm  throughsr"); 


simulation_length=250 .0 ; 

/* 

*/ 

capacity  =  32000; 

/* 

Simulation  Peirameters  */ 

propogation  =  0.7; 

/e 

*/ 

overhead  =  ack.length  =  32.0; 

/e 

*/ 

for  (length  =1024.0;  length  >=  128.0;  length=length/2.0)  { 
for  (logBER  =  -6.0;  logBER  <=  -2.0;  logBER=logBER+ . 1 )  { 

BER  =  pow (10.0, logBER) ;  /♦  ♦/ 

begin_BindoH  =  queue  =  SI  =  i  =1;  /♦  Initialize  variables  ♦/ 

successful.bits  =  buffer  =0.0;  /♦  */ 

while  (i  <=  max.FN)  framestatusCi++]=0; 
m=completed=0: 


top. job  =  newjob(0.0, "send", 0,0.0); 
top_job->previous  =  top.job; 
previous. job  =  top.job; 

job  =  neBjob(simulation_length+l, "send", 0,0.0); 
job->previous  =  previous.job; 
previous. job->following  =  job; 
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/*♦*♦♦♦♦*♦♦***♦♦♦♦♦♦***♦*♦  Begin  Simulating  *********************************/ 
while  (top_job->time  <=  simulation_length){ 


/♦♦*♦♦*♦**♦♦♦♦♦*♦***♦♦♦***♦♦♦♦♦**♦ 


*  Time  to  Transmit  a  Packet 

* 


* 

« 


if  ( ! ( strcmp ( t op_ j  ob->type , " s  end" ) ) ) { 

i  =  begin_windoH ; 
queue  =  SH; 
while  (i  !=  SH)  {. 

if  (framestatus  [i]  ==  3)  •( 
queue  =  i; 
breadc; 

> 

if  (++i  >  max_FN)  i=l: 

> 


/♦  * 

♦  Re-transmit  Packet  in  Queue  ♦ 

*  * 


if  (queue  !=  SN)  { 
tt  =  length/capacity; 

job  =  nowjob(top_job->time+tt,"send",0,0.0) : 
place_job(job) ; 

tr  =  propogation+ ((length+ack_length) /capacity) : 
job  =  newjob(top_job->time+tr, "recv", queue, length) ; 
place_job( job) ; 
framestatus [queue]  =  4; 
framelength [queue]  =  length; 

} 

/♦  ♦ 

*  Transmit  if  Window  not  Full  ♦ 

♦  * 


else  if  (buff er+length  <  max_buffer)  i 
tt  =  length/capacity; 

job  =  newjob(top_job->time+tt, "send” ,0,0.0) ; 
place_job( job) ; 

tr  =  propogation+( (length+ack_length)/capacity) ; 

job  =  newjob(top_job->tiJne+tr, "recv",SN, length) ; 

place_job( job) ; 

buffer=buff er+length ; 

f  r^Lmestatus  [SN]  =  1; 

framelength [SN]  =  length; 

if  (++SN  >  max_FN)  SN  =  1 ; 

> 
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*  IF  Window  is  lull  * 

♦  Try  to  Send  Alter  Receiving  * 

♦  lext  Packet  ♦ 

*  ♦ 

else  ■( 

job  =  nawjob(top_job->lolloiiing->time,''send",0,0.0) ; 
place_job( job) ; 

> 

>  /*  Packet  transmission  complete  */ 


*  ♦ 

*  Time  to  Recieve  a  Packet  4’ 

*  * 

else  il  ( ! (strcmp(top_job->type,"recv'’)))  { 

seed=seed*1103515245+1234S;  /*  Psuedo-rcmdom  */ 

rand=( seed/65636)  */  32768;  /•  generator  ♦/ 

/*  * 

*  HAK  the  received  packet  ♦ 

*  * 

il  (rand/32767  >  powC (1 . 0-BER) ,top_job->length) )  { 

Iramestatus Ctop_job->FN]=3; 

> 

/*  * 

*  ACK  the  received  packet  ♦ 

♦  ♦ 

else  { 

Iramestatus [top_ job->FN] =2 ; 

> 

}  /*  Packet  Received  */ 

♦  ♦ 

♦  Update  Window  Beginning  * 

♦  ♦ 

while  (Iramestatus [begin_window]  ==  2)  •( 

successlul_bits=successlul_bits+lramelength [begin_window] -overhead ; 

Iramestatus  Cbegin_window] =0 ; 

bul 1 er=bul 1 er-lr amelength  Cbegin_window] ; 

il  (++begin_window  >  max_FH)  begin_window=l ; 


♦  * 

♦  Get  next  job  * 

♦  Free  current  job  memory  ♦ 

♦  ♦ 


temp=top_job->lolloBing; 
free (top_ job) ; 
top_job=temp: 

)•  /♦  repeat  lor  next  job  —  *'Hhile(top_job->time<=simtilation_length)"  ♦/ 

/***♦♦♦♦♦♦♦♦♦♦**♦♦*♦**♦•*♦  Simulation  Completed  *******************************/ 

/i^^4‘***rltit**^*****************t**** 

*  * 

*  Record  Throughput  Results  ♦ 

*  * 

ttl  =  lopen("throughsr" ,"a") : 

throughput=successlul_bits/(capacity*simulation_length) ; 
fprintf  (ttf  ,"/i6.5f  y.ll.8f\n", throughput, BER)  ; 
lclose(ttl ) ; 

>  /*  next  BER  --  "lor(logBER=-4.0;logBER<-2.0:logBER=logBER+. 1)"  */ 

)•  /*  lor  (length  =1024.0;  length  >=  128.0;  length=length/2 . 0)  ♦/ 

}  /*  close  main  */ 


52 


*  « 

*  Function  newjob  * 

*  * 
m***************************************************/ 

j  ob_ t  ype  *new j  ob ( t ime , t ype , FH , 1 ength ) 
int  FH ; 

float  time,  length; 

char  type [6] ; 


{ 

job_type  *new: 


if  ( ! (neu=MALLOC ( j  ob_type ) ) )  { 

fprintf (ttf  ,"out  of  the  storage  \n*'): 
f close(ttf ) : 
f close(rf ) ; 
f close(of ) ; 
f close (tf) ; 
exit(l) : 

} 


neB->t ime=t ime ; 
strcpy(neH->type,type) ; 
new->FN=FH ; 
neH->length=length ; 
neB->previous=neB->following=NULL; 
return (new) ; 

> 
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*  * 

*  Fxmction  place_job  * 

*  * 

place_job( job) 
job.type  *job: 

{ 

job_type  ♦current; 
curr ent = t op_ j ob : 

while  (job->time  >=  current->time)  { 
current=current->lollowing ; 

il  ( (current->f ollowing  ==  HULL)  kit  (job->time  >=  current->time) )  { 
curr ent->folloBing= job; 
job->previous=current ; 
job->folloHing=  HULL; 
return ; 

> 

> 

current->previous->folloHing  =  job; 
job->previou8  =  current->previous ; 
curr ent ->prev ious=j  ob ; 
job->f ollouing=current ; 

> 
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/♦  ♦ 

*  Decision  Criteria  for  Level  * 

*  Change  * 

*  * 

’ll  4^  41  ^  if  41  <<  iti  It:  4: 4: 4:  # 4:  <1  It<  #  *  *  4c  / 

if  (m  >=  1)  { 

Ps= (completed-m) /completed ; 
if  (Ps  <  0.87)  { 


length=length/2 . 0 ; 

if  (length  <  128.0)  length=128. 0; 

else  •( 

length  »  length*2.0; 

if  (length  >  1024.0)  length=1024 . 0 ; 

m=completed=0 ; 

} 


Figure  B.l:  Decision  segment  for  single  threshold. 
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/  4<  *  *  4i  *  itc  *  Ik  4i  *  *  4i  If  ^  iti  *  *  lit  *  :•> 


*  * 

*  Time  to  Transmit  a  Packet  ♦ 

♦  * 

If  If  *  Kc  If  If  *  If  If  :k  *  4c  If  %  If  « If  If  If  If  %  If  If  If  If  If  If  #  If  «  He  *  / 

if  (! (strcmp(top_job->type,"send"))){ 

/*  * 

*  Transmit  if  Window  not  Full  * 

★  ♦ 


HcHcifitcHiHcHcitcikik^itiHiHiHc^HiikitcikitcikHcHcHcikHi^HcitcilcHiik/ 

if  (SN-RN  <=  window)  { 


/*  * 

♦  IF  Window  is  full  * 

♦  Go-back-N  ♦ 

♦  * 


4cik4‘4iik4i4iikikikikikik4iikikikikikik4c4!ikikik4cikit'ik4iik4citc/ 

else  { 

SN=RN=1; 

job  =  newjob(top_job->time,"send" ,0 ,0 . 0) ; 
place. job( job) ; 

} 

y  /*  Packet  transmission  complete  ♦/ 


Figure  B.2:  Go-back-N  implementation. 
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1 


> 


/  #  %  41  %  itc  4i  4i  itc  :ti  4c  it>  4i  i|c  :tc  #  4: 4:  it<  If  4^  4^  4<  4c  4c 

4c  4c 

*  Time  to  Recieve  a  Packet  * 

♦  4c 
4c  4c  4c  4c  4c  4c  4c  4c  4>  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  / 

if  ( ! (strcmp(top_job->type,"recv")))  { 
seed=seed*1103515245+12345 ; 
rand=  (seed/65536)  7.  32768; 


/*  * 

*  ACK  the  received  Packet  ♦ 

4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  / 

if  (raiid/32767  <  pow((l  .0-BER)  ,top_job->length))  { 

successful_bits=succ6ssful_bits+top_job->length-overhead; 
if  (SN  ==  1)  SN=0; 
else  SN=1; 

} 

tr* ((length+ack_length) /capacity) +propogat ion; 
job  =  newjob(top_job->time+tr,"recv" ,SN, length) ; 
place_job(job) ; 

} 

temp=top_job->following; 
fr8e(top_job) ; 
top_job=temp; 

>  /4c  repeat  for  next  job  */ 


Figure  B.3:  Stop-and-wait  simulation  block  does  not  require  a  ’’Time  to  Transmit 
a  Packet”  section. 


f 


57 


REFERENCES 


1.  Richard  A.  Yost,  ”011  Nonuniform  Windowing  M-ary  FSK  Data  in  a  DFT 
Based  Detector,”  IEEE  Trans,  on  Communications,  Vol.  COM-28,  No.  12,  pp. 
2014-2019,  December  1980 

2.  Y.  Ishibashi  and  A.  Iwabuchi,  “Throughput  Analysis  of  Adaptive  ARQ  Schemes,” 
Institute  of  Electronics  and  Communications  Engineers  of  Japan,  Vol.  71B,  No. 
6,  pp.  698-707,  June  1988 

3.  S.  Stein,  “Unified  Analysis  of  Certain  Coherent  and  Noncoherent  Binary  Com¬ 
munication  Systems,”  IEEE  Trans,  on  Information  Theory,  Vol.  IT-10,  pp. 
43-51,  January  1964 

4.  I.S.  Gradshteyn  and  I.M.  Ryzhik,  Table  of  Integrals,  Series  and  Products,  pp. 
718  and  961,  Academic  Press,  New  York,  1980 

5.  Harry  L.  Van  Trees,  Detection,  Estimation,  and  Modulation  Theory,  pp.  394- 
396,  John  Wiley  and  Sons  Inc.,  New  York,  1968 

6.  Philip  S.  Yu  and  Shu  Lin,  “An  Efficient  Selective- Repeat  ARQ  Scheme  for  Satel¬ 
lite  Channels  and  Its  Throughput  Analysis,”  IEEE  Trans,  on  Communications. 
Vol.  COM-29,  No.  3,  pp.353-363,  March,  1981 

7.  Kurtis  J.  Guth.  “An  Adaptive  ARQ  Strategy  for  Packet  Switching  Data  Com¬ 
munications  Networks,”  Master’s  thesis.  Naval  Postgraduate  School.  Monterey. 
California,  1989 

8.  David  T.  Tsuda,  “.Adaptive  Go-Back-N:  an  ARQ  Protocol  for  a  Tactical  VS.AT 
Network,”  Master’s  thesis.  Naval  Postgraduate  School,  Monterev.  California, 
1989. 

9.  Dimitri  Bertsekas  and  Robert  Gallager,  Data  Networks,  pp.  50-102,  Prentice- 
Hall,  New  Jersey,  1987 


58 


INITIAL  DISTRIBUTION  LIST 


No.  of  Copies 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5002 

3.  Chairman,  Code  EC  1 

Department  of  Electrical  and 

Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4.  Superintendent,  Naval  Postgraduate  School  5 

Atten:  Professor  Tri  T.  Ha,  Code  EC/Ha 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Superintendent,  Naval  Postgraduate  School  1 

Atten;  Professor  Glen  A.  Myers,  Code  EC/Mv 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  Commandant,  G-PTE-1  2 

U.S.  Coa.st  Guard 

Washington,  D.C.  20593-0001 

7.  Commandant,  G-TPP-2  1 

U.S.  Coast  Guard 

Washington,  D.C.  20593-0001 


59 


8.  U.S.  Department  of  Transportation  Library 
Code  M-493.3,  Room  2200 
400  7th  Street  South-West 
Washington,  D.C.  20590 


60 


